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We present detailed theoretical calculations of two color, time-resolved pump-probe differential 
reflectivity measurements. The experiments modeled were performed on InMnAs/GaSb heterostruc- 
tures and showed pronounced oscillations in the differential reflectivity as well as a time-dependent 
background signal. Previously, we showed that the oscillations resulted from generation of coherent 
acoustic phonon wavepackets in the epilayer and were not associated with the ferromagnetism. Now 
we take into account not only the oscillations, but also the background signal which arises from 
photoexcited carrier effects. The two color pump-probe reflectivity experiments are modeled using 
a Boltzmann equation formalism. We include photogeneration of hot carriers in the InMnAs quan- 
tum well by a pump laser and their subsequent cooling and relaxation by emission of confined LO 
phonons. Recombination of electron-hole pairs via the Schockley-Read carrier trapping mechanism 
is included in a simple relaxation time approximation. The time resolved differential reflectivity 
in the heterostructure is obtained by solving Maxwell's equations and compared with experiment. 
Phase space filling, carrier capture and trapping, band-gap renormalization and induced absorption 
are all shown to influence the spectra. 

PACS numbers: 75.50.Pp, 85.75.-d 



I. INTRODUCTION 

Dilute magnetic semiconductors made of (III,Mn)V 
materials have recently garnered much attention owing to 
the discovery of carrier mediated ferromagnetism. This 
offers the intriguing possibility of their use in the devel- 
opment of semiconductor spintronic devices capable of 
simultaneously performing information processing, data 
storage, and communication functionsi^ 

Time-independent optical studies such as cyclotron 
resonance of III-V DMS materials have provided im- 
portant information on effective masses and electronic 
structure^ exchange parameters j*!* 6 " carrier densities, 7 
and whether carriers are localized or itinerant. 8 

Time-dependent optical studies are even more use- 
ful and can provide information that static magnetiza- 
tion or electrical transport measurements can not. In 
this paper, we report on our theoretical calculations 
and modeling of femtosecond, time-resolved differential 
reflectivity spectroscopy of InMnAs/GaSb heterostruc- 
tures. Femtosecond transient reflectivity spectroscopy 
has proven useful in studying carrier dynamics in semi- 
conductors as well as the generation and propagation of 
coherent phonons in a number of materials. In partic- 
ular, coherent optical phonons have been observed in 
bulk semiconductor s] 9 ! 10 and coherent acoustic phonons 
have been detected in InGaN/GaN-based semiconductor 



heterostructures i 11 : 12 : 13 ; 14 

The experiments we model are two-color pump-probe 
differential reflectivity measurements. In these experi- 
ments, there are several changes to the reflectivity on 
different time-scales. On the fast time scale, there are 
changes to the reflectivity associated with ultrashort car- 
rier lifetimes (~ 2 ps) and multi-level dynamics. 

In addition, pronounced oscillations are observed on a 
longer time scale (~23 ps period) Similar behavior 
was seen in InGaMnAs systems^ Originally these oscil- 
lations were thought to be associated with the ferromag- 
netism in the InMnAs layer since oscillations were not 
observed in samples without Mn. However, we showed 
in a previous work 16 that the oscillations instead resulted 
from selective photoexcitation in the InMnAs layer which 
triggered a coherent phonon wavepacket that propagated 
into the GaSb layer. 

In this paper, we expand upon our previous work and 
show that the large strength of the coherent phonon os- 
cillations results from a fortuitous strong dependence on 
strain of the GaSb dielectric function near the probe en- 
ergy. In addition, we also model the fast time depen- 
dent background signal which arises from the photoex- 
cited carriers. There are three main contributions to this 
transient background signal: 1) the enhanced Drude ab- 
sorption resulting from the photoinduced increase in free 
carriers, 2) the relaxation dynamics associated with the 
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decay of the highly nonequilibrium photoexcited carrier 
distribution and 3) the trapping and subsequent non- 
radiative recombination of photoexcited carriers due to 
the high density of defects in the InMnAs layer. 

We model the experiments by first calculating the de- 
tailed electronic structure in the InMnAs layer. We 
than use a Boltzmann equation formalism to account for 
the photoexcited carrier dynamics. We include band- 
gap renormalization, carrier-phonon scattering, and 
carrier-trapping/recombination through the Shockley- 
Read mechanism. To determine the optical properties, 
we solve Maxwell's equations in the heterostructure. De- 
tails of the experiments and the calculations are given in 
the following sections. 



II. EXPERIMENT 
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FIG. 1: Schematic diagram of the heterostructure consisting 
of a 25 nm Ino.91Mno.09As quantum well and a 820 nm GaSb 
barrier grown on a GaAs substrate. The band gap as a func- 
tion of position in each of the layers is also shown. 



The main sample studied (shown schematically in 
Fig.nj was an InMnAs/GaSb heterostructure, consisting 
of a 25 nm thick magnetic layer with Mn concentration 
0.09, grown on a 820 nm thick GaSb buffer layer on a 
semi-insulating GaAs (001) substrate. Its room temper- 
ature hole density and mobility were 1.1 x 10 19 cm~ 3 
and 323 cm 2 /Vs, respectively, estimated from Hall mea- 
surements. Detailed growth conditions and sample infor- 
mation can be found in Rcf. la 

We performed two-color time-resolved differential re- 
flectivity spectroscopy using femtosecond midinfrared 
pump pulses (2 /xm, ~ 140 fs) and a white light con- 
tinuum probe (0.5 - 1.4 ^m, ~ 340 fs). Experimental 
details are described in Ref. [ljj. 

The source of intense MIR pulses was an optical para- 
metric amplifier (OPA) pumped by a Ti: Sapphire-based 
regenerative amplifier (Model CPA-2010, Clark-MXR, 
Inc., 7300 West Huron River Drive, Dexter, MI 48130). 
At the pump wavelength (2 //m), the photon energy (0.62 
eV) was just above the band gap of InMnAs, so the cre- 
ated carriers had only a small amount of extra kinetic en- 
ergy (~ 0.2 eV at 15 K), minimizing contributions from 
intervalley scattering and intraband relaxation. A white 
light continuum generated by focusing a small fraction 
of the CPA pulses into a sapphire crystal was used as a 
probe, which allowed us to probe a wide energy range 
far above the quasi-Fermi level of the optically excited 
carriers. 



III. THEORY 



A. Bulk InMnAs Bandstructure 



we find it convenient to separate the 8 Bloch basis states 
into upper and lower sets which decouple at the zone 
center. The Bloch basis states for the upper set are 
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which correspond to electron spin up, heavy hole spin 
up, light hole spin down, and split off hole spin down. 
Likewise, the Bloch basis states for the lower set are 



(2a) 
(2b) 
(2c) 
(2d) 



corresponding to electron spin down, heavy hole spin 
down, light hole spin up, and split off hole spin up. 

The effective mass Hamiltonian in bulk zinc blende ma- 
terials in the axial approximation is given explicitly by 21 
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In our pump-probe reflectivity experiments carriers are 
created in the InMnAs layer by pumping below the GaSb 
bandgap. We treat the photogenerated carriers (elec- 
trons and holes) in an 8 band k ■ p effective mass model 
which includes conduction electrons, heavy holes, light 
holes, and split-off holes. Following Pidgeon and Brown 20 



Hn = 



Hi,, Hn 



(3) 



where H uu , H u i, Hi u , and Hu are 4x4 submatrices. The 
effective mass Hamiltonian matrix elements between the 



3 



upper set basis states in Eq. are 
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while the Hamiltonian matrix elements between the lower 
set basis states in Eq. are given by 
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The submatrices coupling upper and lower set basis 
states arc 
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In Eqs. Q -Eg is the bulk band gap, A is the spin orbit 
splitting and k± — k x ±i k y . Following Jain et a&, band 
gap renormalization in InAs shrinks the InAs band gap 
by an amount 



AE g = A N 1/3 



B N 1/4 + C VN 



(5) 



where N is the electron or hole carrier concentration and 
A, B, and C are material parameters which are found 
to be different for electrons and holes. The temperature 
dependence of the band gap is taken into account using 
the empirical Varshni formula 23 . The temperature and 
carrier dependent band gap is given by 



Eg(T,N)=Eg- 
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where a v and f3 v are Varshni parameters which are 
tabulated in Rcf. for a variety of semiconduc- 
tors. The Kane momentum matrix element, V — 
(— ih/mo)(S\p x \X), is related to the optical matrix el- 
ement E p by 2 ^ 
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where k?, = k% + ky 



In the axial approximation 2 ^, we 
have 7 = (272 + 373 )/5 so that the energy bands depend 
only on the magnitude of k|| . Note that at k z — 0, the 
effective mass Hamiltonain in Eq. © is block diagonal 
since L = at k z = 0. 

In Eq. ||SJ| , the parameters 71 and 7 are not the usual 
Luttinger parameters since this is an eight-band model, 
but instead are related to the usual Luttinger parameters 
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This takes into account the additional coupling of the 
valence bands to the conduction band not present in the 
six band Luttinger model for the valence bands. 

The parameter 74 is related to the conduction-band 
effective mass to* through the relation^ 
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The exchange interaction between the Mn ++ d elec- 
trons and the conduction s and valence p electrons is 
treated in the virtual crystal and molecular field approx- 
imation. The resulting Mn exchange Hamiltonian is2£ 



H M n = x N (S z ) 



D a 






-D„ 



(12) 



where x is the Mn concentration, Nq is the number of 
cation sites in the sample, and (S z ) is the average spin 
on a Mn site. The 4x4 submatrix D a is 
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where a and (3 are the s-d and p-d exchange integrals. 

The average Mn spin (S z ) in the ferromagnetic 
Ini-zMn^As quantum well is computed in the mean field 
approximation, i.e. 



(S z ) = -S B S 



JS(S Z 
kT 



(14) 



where Bs{x) is the Brillouin function, S = 5/2 for the 3d 5 
electrons of the the Mn++ ior3, J = 3k B T c /S{S + 1) is 
the ferromagnetic coupling, and T c is the experimentally 
measured Curie temperature. 

The effective mass Hamiltonian for Ini-^Mn^As is 



H = H + H 



Mr. 



(15) 



It is assumed in our calculations that the compensation 
arises from As antisites and hence the effective Mn frac- 
tion, x, in Eq. i|12|) is taken to be equal to the actual Mn 
fraction in the sample. We note that this is supported 
by experimantal evidence showing that InAs grown at 
low temperature (200 °C) is a homogeneous alloy and 
that the magnetization varies linearly with Mn content, 

^ 29.30.31.32 



B. Confined states in the InMnAs quantum well 

In quantum-confined systems such as the InMnAs 
quantum well shown in Fig. ^ we must modify the 
bulk Hamiltonian in Eq. (|15(l . The quantum well breaks 
translational symmetry along the z-direction but not in 
the x-y plane. Since the pump pulse is below the GaSb 
band gap, all photogenerated electrons and holes are 
strongly confined to the well and we assume the con- 
finement potentials are infinite. The wave functions in 
the envelope function approximation are 



(16) 



where A is the cross sectional area of the sample, n is the 
subband index, ky = (fccosfl, fcsin#, 0) is the two di- 
mensional wave vector and \v) = |1) . . . |8) are the Bloch 
basis states defined in Eqs. and (J2J. The complex 
valued envelope functions -Fn^.kn (z) are slowly varying 
in comparison with the Bloch basis states. 

The envelope functions satisfy a set of effective-mass 
Schrodinger equations which, in the axial approximation, 
are 



v' = \ 



^W(ky) F nyM {z) = E n {k) F n ) (17) 



subject to the boundary condition that the envelope func- 
tions vanish at the walls of the quantum well. The opera- 
tors H v v > (k|| ) depend on the wave vector in the x-y plane 
and can be obtained from the matrix elements in Eq. i|15|) 



by making the operator replacement k z —i d/dz in all 
the matrix elements of Eq. J3J). 

In practice, we solve for the envelope functions for a 
given value of ky on an evenly spaced mesh of points, 
Zi, i — 1...N, in the quantum well. Approximating 
the derivative, d/dz by a finite difference formula, the 
Schrodinger equation (|17f) with the rigid wall boundary 
conditions becomes a matrix eigenvalue problem which 
can be solved for the eigenvalues £ n (ky) and the com- 
plex envelope functions -Fn^ky {zi) evaluated at the mesh 
points, Zi. 



C. Boltzmann transport equations 

In the two-color time resolved differential reflectivity 
experiments the pump laser is used to promote electrons 
from the valence to the conduction subbands of the quan- 
tum well. The photoexcited carriers then relax through 
scattering, changing the optical properties of the het- 
erostructure in the process. These processes are often 
simulated using Boltzmann transport equations. 

In this paper, we formulate and solve the Boltzmann 
transport equations using a numerical method similar to 
the one described in Ref. y3 For each subband state 
with energy, E n {k), we have a time dependent distribu- 
tion function, / n (kn,i), which gives the probability, as a 
function of time, of finding an electron in subband n with 
wave vector ky . The Boltzmann equation including pho- 
togeneration of hot electron-hole pairs by the pump, the 
subsequent cooling of these carriers by emission of con- 
fined LO phonons, and the recombination of electron-hole 
pairs by means of carrier trapping is 
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The scattering rate due to scattering by confined LO 
phonons in the quantum well, W k ' n k , , is the rate at which 

electrons in subband n with wave vector ky scatter to 
subband n! with wave vector ky . The last term on the 
right hand side of Eq. I|18l) describes the change in the 
electron distribution function due to the action of the 
pump as well as recombination of electron-hole pairs by 
means of carrier trapping. Thus 
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where the first term on the right hand side is the photo- 
generation rate and the second term is the recombination 
rate due to carrier trapping. 

To simplify the calculations, we use an axial approxi- 
mation in which the distribution functions are replaced 
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by their angular averages in the x-y plane of the quantum 
well. The axial distribution functions are 
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where the limits of integration are from k m — Ah/2 to 
k m + Afc/2. 



Next, we divide fc space into evenly spaced cells of width 
Afc = k max /Nk where Nk is the number of cells. The 
value of k at the midpoint of each cell is denoted k m 
(m = 1 . . . Nk). In each k cell we define the cell averaged 
distribution functions 
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If we assume the distribution functions in the Boltz- 
mann equation depend only on k and vary slowly within 
each fc cell, we can obtain a coupled set of ordinary dif- 
ferential equations for the cell averaged axial distribution 
functions 
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The confined LO phonon scattering rates in the origi- 
nal Boltzmann equation i|18|) depend only on the angle 
6 between ky and kj,, so the cell and axially averaged 
scattering rates appearing in Eq. (|22|l are given by 
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Given initial values of the distribution function, 
fn{k m ,t — > -co), we can solve the system of coupled 
ordinary differential equations (12211 with an adaptive 
Runge-Kutta routine^ 4 In order for the integration to 
be numerically stable, we need the cell averaged scat- 
tering rates in Eq. to satisfy the detailed balance 
condition. We insure this by calculating the downward 
scattering rate and using the detailed balance condition 
to obtain the upward scattering rate. If lf™'^, is the cell 
averaged downward scattering rate between a higher ly- 
ing state E n (k m ) and a lower lying state E n i(k m '), then 
the upward scattering rate, W^,"^, between these two 
states is 
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(24) 



The electrons are initially in thermal equilibrium with 
the lattice and are described by a Fermi-Dirac distribu- 
tion 

1 



fn(kv 
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(25) 



If n and p are the initial electron and hole column 
densities in the quantum well, then the Fermi energy as 
a function of the lattice temperature and initial carrier 
concentrations can be easily found by solving 



(26) 



for Ep using a root finding routined In Eq. the 
valence delta function 8 n>v is defined to be zero if subband 
n is a conduction subband and one if subband n is a 
valence subband. 

To finish specifying the transport problem, we need 
to supply the the cell averaged axial scattering rates ap- 
pearing in Eq. (|22|l as well as the cell averaged photo- 
generation and carrier trapping term. 

D. Photogeneration rates 

The photogeneration rate is computed using Fermi's 
golden rule. The pump is characterized by the fluence 
which is the total flux integrated over time. Assuming 
a narrow spectral width for the pump, centered on the 
pump energy fuv, the fluence is given by 



U (t)- 

n r 



(27) 



where Uo(t) is the pump energy density and n r is the 
refractive index in the InMnAs quantum well. For the 
pump energy density, we assume a Gaussian pulse shape 
with an intensity FWHM of t v so that 
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(28) 



Assuming the Dirac delta function in Fermi's golden 
rule is the only rapidly varying quantity in a k-cell, the 
cell averaged photogeneration rate in Eq. (|19|l is given 
by 
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The cell averaged delta function, denoted A™ n ,(e), is 
dk k 



f (e) 



Ak k 



5 7 {\E n ,{k)-E n {k)\-e) (30) 



where J 7 (x) is a Lorentzian lineshape with FWHW, 7. 

The squared optical matrix element \P n .n'(k)\ 2 be- 
tween subbands n and n' is the angular average 
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where e is the unit complex polarization vector of the 
pump pulse. The optical matrix element between sub- 
band states at fixed kn is 



The electron-hole pair column density is equal to the 
column density of photogenerated electrons which is given 
by 
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where the prime on the summation sign is a reminder 
that the sum over subband index, n, is restricted to con- 
duction subbands. We note in passing that the number of 
photogenerated electrons and holes remain equal in the 
Shockley-Read recombination model. 



F. Confined LO phonon scattering rates 



P„,„'(k||) = — X>|p|«/> / dzF* n ^{z)F n ,, v ,^{z) 

(32) 

where (f|p|i/) is the momentum matrix element be- 
tween the Bloch basis states defined in Eqs. and 
©. Explicit expressions for the matrix elements of 
P = (h/mo) p between the Bloch basis states in terms of 
the Kane momentum matrix element, V, defined in Eq. 
Q can be found in Appendix B of Ref. |jj 



E. Recombination rates due to carrier trapping 

Electron-hole pairs can recombine through the ultra- 
fast trapping of electrons (by Asca antisite defects) and 
holes (by Ga vacancies) by the mid-gap states introduced 
by LT-MBE growthi£ii&. The cell averaged carrier trap- 
ping rate in Eq. (|19l) is treated using a simple relaxation 
time model of the form 
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where /°(fc m ) are the initial thermal equilibrium dis- 
tribution functions defined in Eq. I|25l) and r(t) is 
the Shockley-Read recombination time for electron-hole 
pairs. 

If we assume electron-hole pairs recombine through 
trapping at monovalent flaws in the mid-gap region, the 
Shockley-Read recombination time can be expressed as^ 
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where n and p are the initial electron and hole col- 
umn densities in the quantum well and n e {t) is the col- 
umn density of photogenerated electron-hole pairs. In 
our model, we assume for simplicity that the flaws are 
acceptor-like with t^ <C tq so that ps 0. 



When the sample is excited by the ultrafast pump 
laser, hot carriers are created above the fundamental 
gap. The hot carriers relax back to the band edge and 
reach a quasi-thermal equilibrium through carrier cool- 
ing. The dominant cooling mechanism is absorption and 
emission of confined longitudinal optical (LO) phonons 
in the quantum well. The Frolich Hamiltonian for LO 
phonon scattering in a quantum well of width L is given 
by2£ 
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where the electron-LO phonon coupling constant is given 
by 



ChO = A/47T tuV LO 
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The LO phonon energy in the Einstein model is Hu>lo 
and £0 an d £00 are the static and high frequency dielectric 
constants in the quantum well. The operator ; creates 
a confined LO phonon in the quantum well in the Ith 
LO phonon mode with wavevector q. The vibrational 
amplitude of the ith LO phonon mode is ui{z) and 
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The vibrational amplitude is model dependent, 
slab mode model of Fuchs and Kleiwer— 
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In the 
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Using Fermi's golden rule, the confined LO phonon 
scattering rate due to emission or absorption of a single 
confined LO phonon is 
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The energy slitting A£™£ = - E n >{k') and the 

LO phonon occupation number, No, is given by the Bose- 
Einstein distribution 



Nn = 
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exp(?iw L o/fcsJ 1 ) - 1 
The vibrational amplitude form factor is defined as 
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(41) 



The cell averaged scattering rates, W m ' m , used in the 
Bolztmann equation l|22|l are obtained by substituting 
Eq. (|40() into Eq. (|23|l and performing the integrals. 



G. Generation and propagation of coherent 
acoustic phonons 

The ultrafast photogeneration of electrons and holes 
in the InMnAs quantum well by the pump gives rise 
to coherent longitudinal acoustic (LA) phonons which 
propagate into the sample ^Li2iii2LSi Coherent acoustic 
phonons, as opposed to incoherent phonons, give rise to a 
macroscopic lattice displacement. Since the photogener- 
ated carrier distributions are functions of z, the transient 
lattice displacement U(z,t) due to photogenerated carri- 
ers is independent of x and y and is parallel to z. As 
discussed in Refs. 1371 and flll the coherent phonon lattice 
displacement satisfies a loaded string equation. In the 
presence of a position dependent longitudinal acoustic 
sound velocity, C s (z), we have 



d 2 U(z,t) d 



dt 2 



dU(z,t) 
dz 



where S(z,t) is the driving function, 
acoustic sound velocity is given by 



= S(z,t) (43) 



The longitudinal 



C s (z) = 



Po(z) 



(44) 



where C\\(z) and po(z) are the position dependent elastic 
stiffness constant and mass density. 

The loaded string equation is to be solved subject to 
the initial conditions 

U(z,t = -co) = dU ^ t = -°° ) =0. (45) 

We solve the loaded string equation numerically by fi- 
nite differencing Eq. I|43|) inside a computational box 



whose left edge, zl, is the semiconductor-air interface 
and whose right edge, zr, lies inside the GaAs substrate 
(see Fig. QJ. At zr we impose absorbing boundary con- 
ditions and at there are no perpendicular forces at 
the semiconductor-air interface. Thus we solve the ini- 
tial value problem subject to the left and right boundary 
conditions 



dU(z L ,t) 



dz 







and 



dU(z R ,t) 
dz 



Cs(zr) 



dU{z Rl t) 
dt 



0. 



(46a) 



(46b) 



Starting with the second quantized Hamiltonian for 
the electron-phonon interaction, a microscopic expression 
for the driving function was derived in Ref. |37] using the 
density matrix formalism (see also the erratum in Ref. |3^ 
as well as the review article in Ref. ITlTl . In zinc-blende 
materials such as InAs the electron-phonon interaction 
is due to deformation potential coupling. Under typical 
experimental conditions the microscopic expression for 
the driving function can be simplified to^i 



S(z,t) = - 
Po 



dp e (z,t) dp h (z,t) 



dz 



dz 



(47) 



where po is the mass density, a c and a v are the defor- 
mation potentials for conduction and valence bands and 
p e (z,t) and ph(z,t) are the photogenerated electron and 
hole carrier densities. We note that this last equation was 
derived independently in the elastic continuum limit by 
Chigarev et al22i and Chern et al»ii The driving function 
satisfies the sum rule 



dz S(z,t) = 0. 



(48) 



as shown in Refs. |33 andllll 

The photogenerated electron and hole densities can be 
obtained from the envelope and distribution functions 
described in the previous sections as 



n,!Aki| 



PeM) = 7 E {fn(k,t)-f°(k))\F n ^(z)\ 2 5 n , c 

(49a) 

Ph (z,t) = j (fn(k)-fn(k,t))\F niViH (z)\ 2 S n , v . 
n,^,ki| 

(49b) 



and 
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FIG. 2: Background model dielectric function at T = K as 
a function of photon energy used in computing the dielectric 
function in the InMnAs quantum well. 



The initial Fermi-Dirac distribution functions /„(&) are 
defined in Eq. (|25|) and S nyC and S n>v are conduction and 
valence delta functions that select out conduction and 
valence subbands, respectively. 

The propagating coherent phonon displacement field 
U(z,t) gives rise to a propagating strain tensor with com- 
ponents 



£zz{z,t) 



dU(z,t) 
dz 



(50a) 



function in the entire structure, i.e. the InMnAs well, 
GaSb barrier, and GaAs substrate, over the probe energy 
range which extends up to 2 eV. We will denote this 
dielectric function as 

e{fuo, z, t) = e±(huj, z,t) + i S2(huj, z, t) (51) 

where fku is the probe energy. Once the dielectric func- 
tion is found, we can solve Maxwell's equations for the 
time dependent probe reflection coefficient, R(huj, t), and 
the transmission coefficient, T(hu>,t), using the transfer 
matrix method described in detail in Ref. |4£j. 

There are several processes which contribute to the di- 
electric function in the InMnAs quantum well. The first 
of these is a Drude term due to free carriers in the quan- 
tum well which gives a real contribution to the dielectric 
function of 



(M 2 



(52) 



where u> p is the plasma frequency. The Drude contribu- 
tion to the dielectric function in Eq. (|52|l is uniform in 
the quantum well and vanishes everywhere else. In the 
random phase approximation (RPA) , the time dependent 
plasma frequency is given byai 



and 



Aire 2 



L A 



fn(k,t)-S 

1 1 . 

m*(fc) 



(53) 



C 12 (z) 



Cn(z) + C 12 (z) 



£zz{z,t) 



(50b) 

where Cn(z) and C\i{z) are elastic stiffness constants. 
This propagating strain field alters the optical properties 
of the sample which can be detected by the probe. 



where L is the quantum well width and A is the cross 
sectional area of the sample. The effective mass is 



m*(fc) 



1_ d 2 E n (k) 
h 2 dk 2 ' 



(54) 



H. Transient probe response 

To compute the time dependent probe transmission 
and reflection coefficients we need to model the dielectric 



A second contribution to the dielectric function in the 
quantum well is due to dipole transitions between quan- 
tum confined carrier states. Using Fermi's golden rule, 
we obtain22i2£ 



J 



Si{hjJ,Z,t) QW = "^ \ P n,n'{k\\)\' 



/n'(M)-/»(M) 



L A ^ 

n.n' ,feii 



AE 7hn ,(k) (AE 7hn ,(k) + Hlu) (A£„, n /(fc) - Huj) 



(55a) 



and 



e 2 (hj,z,t) QW 



47r 2 e 2 1 

L (M 2 A 



J2 l P »,n'(fc||)| 2 (/n'(M)- /«(M)) 5{AE n ^{k)-hw) 



(55b) 
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where AE n<n r(k) = E n (k) — E n /(k) are the transition 
energies, including band gap renormalization corrections 
due to photogenerated carriers, and P n ,n'(k\\) are the op- 
tical matrix elements. The contributions to the dielectric 
function in Eq. (|55|l are for zero linewidth. For a finite 
FWHM linewidth T, we make the replacements^ 



A£ n>Ti /(fc) - fiw (AE n , n ,(k) - hLo f + {r/2y 



and 



S(x) 



(r/2) 



7T x 2 + {T/2y 



(56) 



(57) 



There is also a background dielectric function, e b in the 
quantum well due to all the higher lying electronic tran- 
sitions whose real and imaginary parts we shall denote 
e%t, and £26- For simplicity, we treat these contributions 
to the dielectric function using Adachi's model dielectric 
function for bulk InAs with contributions from the Eq 
and Eg + Ag critical points removed424i These corre- 
spond to contributions from the confined quantum well 
electronic states and are already included in Eq. 1)55(1. 
The background dielectric function in InAs as a function 
of photon energy is shown in Figure |21 at T = K and in 
the absence of strain. 

Following Thomsen et. al4^& we assume that the 
dielectric function changes with strain only because of 
strain induced variations in the energy gaps associated 
with each transition. The propagating coherent phonon 
strain tensor alters the optical properties of the structure 
through the deformation potential interaction. In our ex- 
periments, the probe photon energy can go as high as the 
GaSb Ei transition region. Ab initio density functional 
calculations of the deformation potentials for the E\ tran- 
sitions in a number of semiconductors^ have shown that 
the deformation potentials associated with the Eq and 
Ei features are equal to within 20%. So to a first ap- 
proximation the effect of temperature and strain on e b is 
to introduce a rigid shift in the dielectric function such 
that 

e b (fuv, z, t) = e b (hw - AE g (T) - a cv (e xx + e yy + e zz )). 

(58) 

Here, a cv = a c — a v is the interband deformation po- 
tential, e xx , e yy , and e zz are the coherent phonon strain 
tensor components defined in Eq. JSUJ), and AE g (T) — 
E g {T) — Eg is the band gap shift due to temperature 
variations with E g (T) being the temperature dependent 
band gap defined by the Varshni expression in Eq. JSJ . 

The total dielectric function in the quantum well is 
obtained by adding the Drude, Quantum well, and back- 
ground contributions in Eqs.JS^J), l|55|) and l|58|) . i.e. 

e(huj, z,t) = £(huj, z,£)d + e(huj, z, t)Qw + £b(hu, z, t) 

(59) 

We use Adachi-type model dielectric functions for the 
GaSb barrier^24^ and the GaAs substrate^ Fig. [3]shows 
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FIG. 3: Model dielectric functions at T = K as a function 
of photon energy for (a) GaSb and (b) GaAs used in calcu- 
lating the dielectric functions in the GaSb barrier and GaAs 
substrate. 



the real and imaginary parts of the model dielectric func- 
tion for bulk GaSb and GaAs at T — K in the absence 
of strain. Temperature and strain effects in the GaSb 
barrier and GaAs substrate are included using the same 
rigid shift model as defined in Eq. I|58|) . Note that the 
dielectric functions in GaSb and GaAs are modulated by 
the coherent phonon strain field as it propagates through 
the structure. 



IV. RESULTS AND DISCUSSION 

A. Quantum well electronic states 

The 25 nm thick ferromagnetic Ino.91Mno.09As quan- 
tum well is shown schematically in Fig.f^ The quantum 
well is p-type with a free hole density estimated to be 
p « 10 19 cm -3 and a Curie temperature of T c = 55 K. 
The computed bandstructure of the 25 nm ferromagnetic 
Ino.91Mno.09As quantum well in the axial approximation, 
assuming infinite barriers, is shown in Fig. 0] as a func- 
tion of fcy for temperatures above and below the Curie 
temperature. The Fermi energies, Ef, corresponding to 



a free hole density of p — 10 cm~ J arc indicated by 
short horizontal lines in Fig. 0] The bandstructure well 
above the Curie temperature (T = 300 K) is shown in 
Fig. 01a). Far above the Curie temperature, the aver- 
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FIG. 4: Bandstructure of a 25 nm ferromagnetic Ini-zMrizAs 
quantum well with infinite barriers and x = 9 % at (a) T = 
300 K and (b) T = 20 K. The Curie temperature is taken to 
be T c — 55 K. The Fermi energies, Ef, for a hole density of 
p = 10 19 cm -3 are indicated by the horizontal lines. 



age Mn spin (S z ) vanishes, the sample is nonmagnetic, 
and the computed subbands are doubly degenerate. Be- 
low the Curie temperature, the quantum well becomes 
ferromagnetic with a nonvanishing (S z ) and the doubly 
degenerate subbands in Fig.^Ja) become spin-split as can 
be seen in Fig. |H[b). 



B. Coherent phonon generation and propagation 



FIG. 5: Coherent phonon strain field as a function of position 
for delay times ranging from to 140 ps. The curves have been 
offset for clarity. 

quantum well at long times. The fact that the steady- 
state strain is localized in the well follows directly from 
the sum rule 1(4*51) . 

In addition to the localized strain, a transient strain 
pulse propagates into the GaSb barrier at the longitu- 
dinal acoustic sound speed. Two transient strain pulses 
are generated in the well, one propagating to the left and 
the other to the right. The leftward propagating pulse is 
totally reflected off the semiconductor-air interface and 
trails the rightward propagating pulse as it propagates 
into the GaSb barrier. 



The coherent phonon lattice displacement, U(z,t), ob- 
tained from the loaded string equation 143(1 gives rise to 
strain tensor components e xx (z 1 t), e yy (z,t), and e zz (z,t) 
as defined in Eq. I|50|) . The coherent phonon strain ten- 
sor component e zz [z,t) is shown in Fig. [S] as a function 
of position in the sample for several equally spaced delay 
times ranging from to 140 ps. 

Following photogeneration of carriers by the pump, a 
localized strain appears in the quantum well as can be 
seen in Fig. |3J This is due to near steady-state load- 
ing by the driving function at long times. Assuming the 
driving function, S{z,t), is approximately time indepen- 
dent at long times, the loaded string equation l(43|) can be 
integrated once in the steady-state limit. The resulting 
steady-state strain is 

, s dU(z) f z , , S(z') , s 

where C s is the longitudinal acoustic sound speed in the 
InMnAs quantum well and S(z') is the approximately 
time independent driving function left behind in the 



C. Differential reflectivity 

The oscillation observed in the differential reflectiv- 
ity can be attributed to propagation of the strain pulse 
through the structure. During most of the experiment, 
the travelling strain pulse is in the GaSb barrier shown 
schematically in Fig. ^ 

The propagating strain tensor shown in Fig. [S] alters 
the dielectric function in the structure. The change in 
the complex dielectric function due to coherent phonon 
wavepackets is given by 

Ae(fa,z,t)= d£i ^' z) e zz {z,t), (61) 
de zz 

where the total derivative with respect to strain, 
de{hijj, z) I 'de zz , is piecewise constant in z having differ- 
ent values in the InMnAs well, the GaSb barrier, and the 
GaAs substrate. 

The total derivative with respect to strain in Eq. is 
obtained from the Adachi-type model dielectric functions 



11 



CO 
CO 

"a 



200 



100 



-100 



-200 



GaSb / 


Re 




Im 







12 3 4 

Photon Energy (eV) 



FIG. 6: Derivative of the complex GaSb dielectric function 
with respect to strain as a function of the probe photon en- 
ergy. The solid line is the real part and the dashed line is the 
imaginary part. 



by differentiating Eq. (|58|l with respect to e zz taking care 
to eliminate e xx and e yy in favor of e zz using Eq. (|50|) . In 
the case of GaSb, the real and imaginary parts of de/de zz 
are plotted as a function of the probe photon energy in 
Fig. [fj] As seen in Fig. [fj] the perturbation in the dielec- 
tric function per unit strain, i.e. the dielectric function 
strain sensitivity, is a function of the probe energy. In 
particular, there is a giant peak in the dielectric function 
strain sensitivity at the GaSb E\ transition near 2.0 eV. 
Consequently, the best probe wavelength for observing 
coherent phonon differential reflectivity oscillations is in 
the region around the E\ transition. 

The differential reflectivity at a given probe delay 
time is obtained from the dielectric function by solv- 
ing Maxwell's equations in the entire structure using 
the transfer matrix formalism as described earlier. If 
we use the complete time- and space-dependent dielec- 
tric function defined in Eq. I|59|l including the quantum 
well, Drude, and background contributions, we can com- 
pute the total differential reflectivity containing the co- 
herent phonon oscillation and transient carrier relaxation 
effects. 

Our theoretical results are compared with experiment 
in Fig. where we plot the experimental and theoretical 
differential reflectivity spectra for probe wavelengths of 
650 and 775 nm. In both cases there is an initial sharp 
drop in the differential reflectivity which we attribute to 
free carrier Drude absorption by the hot carriers created 
by the pump. 

The photogenerated hot carriers relax back to quasi- 
cquilibrium distributions at their respective band edges 
through emission of confined LO phonons. The relax- 
ation of photogenerated carriers by LO phonons alters 
the quantum well dielectric function in Eq. (|55|l through 
changes in the time-dependent distribution functions. 
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FIG. 7: Theoretical and experimental differential reflectivity 
for probe wavelengths of (a) 650 and (b) 775 nm due to vari- 
ations in the time- and space-dependent dielectric functions. 



This carrier cooling by LO phonon emission results in 
the subsequent rise in the differential reflectivity traces 
seen if Fig. [7| 

In addition to carrier cooling by LO phonon emission, 
electron-hole pairs recombine through trapping at mid- 
gap defects with To « 200 ps. This gives rise to the slow 
decay in the differential reflectivity at long times seen 
in Fig. [7| At short times, electron-hole pair recombina- 
tion is enhanced since the Shockley-Read recombination 
time, 7"(i), is a monotonically decreasing function of the 
photogenerated electron-hole pair density. 

For delay times of less than 20 ps, our theory doesn't 
agree very well with the experiment. For delay times 
greater than 20 ps, however, the theory reproduces the 
experimental results surprisingly well. In particular, the 
period and amplitude of the reflectivity oscillations in 
relation to the height of the plateau as well as the decay 
of the reflectivity oscillations with delay time are in good 
agreement with experiment. 

The oscillations in differential reflectivity seen in Fig.[7| 
are due to changes in the background dielectric function 
induced by the propagating coherent strain pulse seen in 
Fig. [SJ If we compute the probe differential reflectivity 
neglecting the quantum well and Drude contributions to 
the total dielectric function in Eq. I|59|l and retain only 
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FIG. 8: Experimental (a) and Theoretical (b) coherent 
phonon differential reflectivity oscillations for probe wave- 
lengths of 650, 775 and 850 nm. 
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the background contribution, we get the coherent acous- 
tic phonon differential reflectivity oscillation absent the 
transient relaxation signal. 

In Fig. [3] the computed coherent phonon differential 
reflectivity oscillations are shown as a function of time 
delay for probe wavelengths of 650, 775 and 850 nm, cor- 
responding to photon energies of 1.9, 1.6, and 1.46 eV re- 
spectively. The theoretical differential reflectivity curves 
in Fig. El (b) agree well with the experimentally measured 
differential reflectivity seen in Fig. |H1 (a) after subtraction 
of the transient background signal. As we go from 650 
to 850 nm, the differential reflectivity oscillation period 
becomes longer. 

The reflectivity oscillations can be qualitatively under- 
stood as follows. The propagating strain pulse in Fig. 
gives rise to a perturbation in the GaSb dielectric func- 
tion which propagates at the acoustic sound speed. The 
sample thus acts as a Fabry-Perot interferometer and a 
simple geometrical optics argument shows that the period 
for the reflectivity oscillations due to the propagating co- 
herent acoustic phonon wavepacket is approximately^ 



T = 



A 



2 C s n(A) 



(62) 



where A = 2ir c/ui is the probe wavelength, C s is the LA 
sound speed in the GaSb barrier and n(A) is the wave- 
length dependent refractive index. The refractive index 
can be obtained from the GaSb model dielectric function 
in Fig. (a) aa£i 



i(A) = \j\ (ei(A) + ^"i(A) 2 + £2(A) 2 ) 



(63) 



In Fig.[5]we have plotted experimentally measured coher- 
ent phonon differential reflectivity oscillation periods as 
a function of probe wavelength as solid circles. The solid 
line shows the oscillation period vs. probe wavelength es- 
timated using Eq. I|62|l . The excellent agreement between 



FIG. 9: Coherent phonon differential reflectivity oscillation 
period vs. probe wavelength. The experimental data are 
shown as solid circles, the dashed line shows the oscillation 
period calculated from the theory described in the text, and 
the solid line shows the oscillation period estimated from Eq. 



theory and experiment is compelling evidence that the re- 
flectivity oscillations seen in the experiments are induced 
by propagating coherent acoustic phonons in the GaSb 
barrier. 

In going from a probe wavelength of 650 to 850 nm 
in Fig. HO we note that the initial amplitude of the dif- 
ferential reflectivity oscillation decreases with increasing 
probe wavelength. At the same time these oscillations 
become more weakly damped. The reason for the re- 
duction in amplitude of the oscillations can be found in 
Fig. ITUl where we plot de/de zz as a function of probe 
wavelength. As the probe wavelength increases (and the 
photon energy decreases), the strength of the perturba- 
tion of the dielectric function due to the propagating co- 
herent phonon strain defined in Eq. I|61|l decreases. This 
accounts for the observed reduction in the initial ampli- 
tude of the differential reflectivity oscillations as we go 
to higher wavelengths. The increased damping of the 
differential reflectivity oscillations with decreasing probe 
wavelength is simply due to the fact that the absorp- 
tion coefficient in GaSb is rapidly decreasing with wave- 
length in this wavelength range as can be inferred from 
the imaginary part of the GaSb dielectric function plot- 
ted in Fig. (a). 



V. CONCLUSIONS 

In summary, we have performed calculations and mod- 
elled time-dependent two-color differential reflectivity ex- 



13 



Probe Energy (eV) 




500 600 700 800 900 

Probe Wavelength (nm) 



FIG. 10: Derivative of the complex GaSb dielectric function 
with respect to strain as a function of the probe wavelength 
in the experimentally relevant wavelength range. The solid 
line is the real part and the dashed line is the imaginary part. 



periments on a ferromagnetic InMnAs/GaSb heterostruc- 
ture. We have observed large amplitude reflectivity oscil- 
lations resulting from the generation of coherent acoustic 
phonon wavepackets in the InMnAs quantum well and 
their subsequent propagation into the GaSb layer. The 
propagation of these coherent, localized strain pulses into 
the GaSb buffer results in a position- and frequency- 



dependent dielectric function. 

To take into account the time dependent background 
differential reflectivity, we modeled the two color pump- 
probe reflectivity experiments in a Boltzmann equation 
formalism. Electronic structure was calculated using k ■ p 
theory in a confined InMnAs layer. We included 1) pho- 
togeneration of hot carriers in the InMnAs quantum well 
by a pump laser and 2) their subsequent cooling by emis- 
sion of confined LO phonons. Recombination of electron- 
hole pairs via the Schockley-Read carrier trapping mech- 
anism was also included in a simple relaxation time ap- 
proximation. 

Our results agree remarkably well with the experimen- 
tal coherent phonon oscillations and reasonably well with 
the time-dependent background signal and capture the 
major qualitative trends of the data. We identify three 
key effects which contribute to the backgrounds signal: 
1) the enhanced Drude absorption resulting from the in- 
crease in carriers from the laser photoexcitation, (nega- 
tive AR/R), 2)the relaxation dynamics associated with 
the decay of the highly nonequilibrium photoexcited car- 
rier distribution (positive AR/R) and 3) the trapping 
and then non-radiative recombination of the photoex- 
cited carriers resulting from the high density of defects 
in the InMnAs layer (positive AR/R). 
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